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^^ \ Abstract. Properties of Lyman Break Galaxies at z ~ 5 are presented and 

^^ , are compared with those at z ~ 2 — 3 LBGs to see evolution of LBGs from z ^ b 

CN ' to z ~ 3 — 2. The LBGs at z ^^ 5 are tend to have smaller stellar masses, but 

they may show large star formation rates and specific star formation rates. 
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Introduction 



We have been studying properties of Lyman Break Galaxies (LBGs) at z ~ 5, 
~ 1 Gyr after the Big Bang and before z ~ 3. In this contribution, we present 
^ ' our recent results obtained for 2: ~ 5 LBGs and possible evolution from 2: ~ 5 

C3 . LBGs to z ~ 2 — 3 LBGs. Our sample at z ^ b was made in a field including 

the Hubble Deep Field North (HDF-N) by using the Subaru SCam through V-, 
Ic-, and z'-filters. Since in and around the HDF-N field many spectroscopic 
^ '_ data were available, the color criteria to select the z ~ 5 LBGs can be reliably 

examined. A sample of z ~ 5 LBGs was chosen from galaxies brighter than 
z' = 26.5 mag. A number of the sample is ~ 600 in an area of ~ 500 arcmin^. 



f— ^ ■ 2. Luminosity dependent evolution? 

00 ■ 

^D ■ UV luminosity function is presented by Iwata et al. (2003, 2007, 2008 this vol- 

ume). In their results, from z ~ 5 to ~ 3, a number of faint LBGs (fainter than 

, , ■ L*) is increasing with redshift, while in the brighter part the number does not 

r> I change significantly. Although there is a controversy on this result, a luminosity 

C^ ■ dependent evolution of UV luminosity function is indicated. 

Spectroscopic follow-up observations of a part of our LBG sample were made 
by Ando et al. (2004). They found that UV luminous LBGs (brighter than L*) 
do not show large equivalent width of Lya emission line. Ando et al. (2006) 
showed this trend together with other LBGs as well as Lya emitters (LAEs). 
This may be due to the dusty environment or the presence of more massive HI 
gas in the luminous LBGs. 
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3. Stellar masses and star formation rates in z ~ 5 LBGs 

Using public IRAC data in the GOODS-N region and those obtained in the 
flanking fields, we constructed a z ~ 5 LBG sample with IRAC photometry. 
From the sample of the LBGs mentioned above, we selected the LBGs those 
are free from the contamination by neighboring objects by eye inspection. The 
sample size is ~ 170; about half of them are detected both in IRAC channels 
1 and 2. We derived the Spectral Energy Distributions (SEDs) for the sample, 
and made SED fitting. The spectral model adopted for the SED fitting we use 
is rather traditional, because we intend to compare with the results obtained 
in the previous studies at z ~ 2 — 3; stellar synthesis code by Bruzual and 
Chariot (2003) with 0.2 Z©, the Salpeter IMF (0.1 - lOOM©) under constant 
star formation history, and Calzetti extinction curve (Calzetti et al. 2000). We 
include the Ha emission line calculated from the model intrinsic star formation 
rate and extinction. The SED fitting method is the essentially the same as that 
by Sawicki and Yee (1998). In the fitting, we fixed the redshift to be 4.8. 

The resulting stellar masses, star formation ages, color excess, and star 
formation rates are presented in Figure [1] (left) overlayed with those obtained 
for z ~ 3 LBGs (Shapley et al. 2001) and for z ~ 2 LBGs (Shapley et al. 2005). 
The rest-frame UV luminosity and optical luminosity of these samples are mostly 
overlapped each other. The stellar masses of the z ~ 5 LBGs are significantly 
smaller than those at z ~ 2 — 3, and ages are smaller on average. These suggest 
the premature nature of the z ~ 5 LBGs. Stellar mass function and stellar mass 
density at z ~ 5 are presented by Yabe et al. (2008, this volume). 



4. LBGs with large SFRs? 

As seen from Figure [T] (left), the star formation rates (SFRs) are systematically 
larger than those in z ~ 2 — 3 LBGs. Figure [1] (right) shows SFRs against their 
stellar masses for z ~ 5 LBGs (filled circles). The SFRs are very high amounting 
to several hundreds Mq yr~^. These values are significantly larger than those 
seen in z ~ 2 — 3 LBGs shown as open symbols in Figure [1] (right). They 
are also larger compared with the results for z ~ 0.1 — 2 star forming galaxies 
(e.g., Daddi et al. 2007). The specific star formation rates tend to decrease 
with decreasing redshift. It is suggested that the z ~ 5 LBGs are in active star 
forming phase and have dusty environment. In fact, the color excesses obtained 
for the z ~ 5 LBGs tend to be larger than those obtained for z ~ 2 — 3 LBGs 
(Figure [T][left)). We may be witnessing an emergence of a population of active 
star forming galaxies at z '^ 5. 

If the z ~ 5 LBGs keep the high SFRs, they should have stellar masses 
more than lO^^M© at z ~ 3. The distributions of the expected stellar masses at 
z ~ 2 — 3 are shown as dotted lines in Figure mieft), suggesting that the z ~ 5 
LBGs are not the progenitor of the z -^ 2 — 3 LBGs or the large SFRs should 
decline. If the SFRs decline exponentially with a time scale of 10^ yr, SFRs 
should be very much small at z ~ 3. Thus the star formation may be episodic if 
the z ~ 5 LBGs are the progenitor of z ~ 3 LBGs. Or they might be progenitor 
of passive massive galaxies at z = 2 — 3. 
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Figure 1. Left: Stellar mass (upper left), age (upper right), color excess 
(lower left), and star formation rate (lower right) of the SED fitting are shown 
with solid lines. Also shown are those for z ~ 2 and z ~ 3 LBGs (Shapley et 
al. 2005; Shapley et al. 2001). Right: Star formation rate vs stellar mass 
for z ^ 5 LBGs (filled circles). Also shown are those for z ~ 2 — 3 LBGs 
(open symbols). Solid lines show regression lines for star forming galaxies at 
the redshifts shown taken from Daddi et al. (2007). 



However, as known well, there may be an age-extinction degeneracy in the 
SED fitting. Because we do not have NIR photometric data, our SED fitting 
may suffer from the degeneracy. To examine this, we made SED fitting for 
LBGs at z ~ 5 sampled by Stark et al. (2007) in which J and K photometry 
are available. We made SED fitting with and without J and K data. The 
resulting output parameters are shown in Figure [2] (left). The stellar masses are 
robust as seen from the figure. The SFRs are also rather robust (except for a 
few extraordinary cases), though the error bars are large. We also estimated 
the color excess using the method by Meurer et al. (1995), and found many 
of the objects have the color excesses consistent with the SED fitting. We also 
made various tests to see whether the large SFRs are due to some artificial 
causes. Figure [2] (right) shows output parameters of the SED fitting when we 
take redshift as a free parameter (upper left panel), other star formation history 
(upper right panel), metallicity (lower left panel), and other extinction curves 
(lower right panel). As seen from the panels, except for the adoption of other 
extinction curves, majority of estimated SFRs agree with those obtained in the 
fiducial model. 

No LBGs with large SFRs coincide with SCUBA sources (Pope et al. 2005), 
X-ray sources (Alexander et al. 2003), and VLA sources (Richards 2000). The 
estimated upper limits on SFRs from these surveys are ~ 15OOM0 yr~^. It is 
worth noting that many of these LBGs with large SFRs do not reside in the 
SCUBA/Chandra surveys. Further observations are necessary to definitively 
conclude the presence of such population. 
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Figure 2. Uncertainty of SED fitting. Left: stellar mass, age, color excess, 
and star formation rate obtained from the SEDs by Stark et al. (2007) with 
(abscissas) and without (ordinates) J and K data. Right: Star formation 
rates using fiducial model (abscissas) and other models (ordinates). Cases 
by adopting free parameter redshift, star formation history of exponentially 
decaying model (triangles) and two-component model (squares), metallicity 
of 1.0 Zq (circles) and 0.1 Zq (triangles), and extinction curves of the Milky 
Way (squares), the LMC (triangles), and the SMC (circles) are shown. 
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